CLAIMS 



What is claimed is: 

1. A yaw stability system for a vehicle having a plurality of wheels each with a 
torque control element, said yaw stability system comprising: 
a yaw rate sensor measuring a vehicle yaw rate; 

a plurality of braking devices each operably associated with one of the torque control 
elements and configured to exert a braking torque on the control element in response to a 
control command; and 

a control unit communicating with said yaw rate sensor and configured to 
identify a desired yaw rate, 

determine a yaw rate tracking error based on the difference between the desired 
yaw rate and the vehicle yaw rate, 

determine a yaw condition of the vehicle based on the vehicle yaw rate, 
determine a control yaw moment to minimize the yaw rate tracking error, 
select one or more of said plurality of braking devices based on the yaw 

condition, and 

communicate a control command to the one or more selected braking devices to 
induce said control yaw moment. 



2. The yaw stability system of claim 1 wherein said control unit determines the 
control yaw moment using a sliding mode control law based on a lumped mass vehicle model. 
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3. The yaw stability system of claim 2 wherein said control unit determines the 
control yaw moment (Mz) based on the following equation 

M z =IJ des ~[a[(C FL +C fR )a f eo*8 + (ri FL +ri FR )F zF sin S)-b(C RR + C R£ )a R + 

(c*C pL -d* C fR )a f sin S - c(^F f cos S + tj^F^ ) + rffo^ cos * + ^ )] + I 22 ?]SAT(^^) 

4. The yaw stability system of claim 1 wherein said braking devices are 
electromagnetic retarders and wherein said control command is a current command. 

5. The yaw stability system of claim 4 wherein said electromagnetic retarders are 
eddy current machines. 

6. The yaw stability system of claim 5 wherein the control unit estimates braking 
device saturation torque (T est ) based on a quadratic function of rotor speed and excitation 
current. 

7. The yaw stability system of claim 6 wherein the quadratic function is: 

T est -a^f^ri^f 2 (^yi 2 

8. The yaw stability system of claim 7 wherein the coefficient functions /ofco), 
/i(co),and /2(a)) are estimated from steady state test data performed for various rotor speeds. 
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9. The yaw stability system of claim 8 wherein the coefficient functions /o (cg), 
/i(o), and /2(a) are defined by 

f i (o)) = a i0 +a n co + a i2 a) 2 

and wherein the parameters o$j are estimated through a least square fit based on the steady state 
test data. 

10. The yaw stability system of claim 9 wherein the coefficient functions /o(co), 
/i(co), and f 2 ((ri) are estimated by recalculating the coefficient functions fo((o\ /2(a)), and /2(G)) 
for each rotor speed based on the estimates of parameters Oy and the parameters ojj are then 
estimated based on the recalculated coefficient functions /o(co), /i(o), and f 2 ((o) through a least 
square fit based on the steady state test data. 

11. The yaw stability system of claim 4 wherein said control unit is an open loop 
controller providing a current optimal torque without a current feedback signal from the 
electromagnetic retarder. 

12. The yaw stability system of claim 4 wherein the control unit is a parametric 
model control unit. 
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13. A method for controlling yaw in a vehicle having front left, front right, rear left, 
and rear right wheels and a plurality of braking devices each associated with one of the wheels, 
said method comprising: 

determining a vehicle yaw rate; 
determining a desired yaw rate; 

calculating a yaw rate error based on the difference between the desired yaw rate and 
the vehicle yaw rate; 

determining a control yaw moment using a sliding mode control law based on a lumped 
mass vehicle model; 

selecting one of the braking devices based on a vehicle yaw condition; 

determining a control command for the selected braking device based on the control 
yaw moment; and 

communicating the control command to the one or more selected braking devices. 

14. The method of claim 13 wherein the step of determining a control command 
further includes determining a required torque for the selected braking device, said required 
torque being the torque required from the selected braking device to induce the control yaw 
moment. 

15. The method of claim 14 wherein the plurality of braking devices are eddy 
current machines, wherein the control command is a current command, and wherein the step of 
determining the current command further includes determining a saturation torque for the 
selected braking device based on a quadratic function of control element speed and excitation 
current. 
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16. The method of claim 15 wherein the step of determining the current command 
further includes comparing the saturation torque for the selected braking device to the required 
torque. 

17. The method of claim 16 wherein the step of determining the current command 
further includes determining a command current for the selected braking device if the required 
torque is less than the saturation torque. 

18. The method of claim 16 wherein, if the required torque is greater than the 
saturation torque, the step of communicating the current command farther includes sending a 
saturation current command to the selected braking device, selecting a second braking device, 
and sending a second current command to the second braking device to cause the second 
braking device to exert a torque equal to the difference between the control yaw moment and 
the saturation torque, and wherein the first and second selected braking devices are on the same 
lateral side of the vehicle. 

19. The method of claim 15 wherein the step of determining the saturation torque 
(Test) is based on the following equation: 

and wherein the method further includes estimating the coefficient functions /o(co), /i(co), and 
/ 2 (co) from steady state test data performed for various rotor speeds. 
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20. The method of claim 19 wherein the coefficient functions /o(co), /i(co), and 
/2(a)) are defined by 

f i (a)) = a i0 +a n a) + a i2 6) 2 

and wherein the step of estimating the coefficient functions /o(co), /1(a)), and /2(G)) includes 
estimating parameters ay through a least square fit based on the steady state test data, includes 
recalculating the coefficient functions /o(co), /2(w), and /2(a)) for each rotor speed and based on 
the estimates of parameters ajj, and re-estimating the parameters QSj based on the recalculated 
coefficient functions /o(co), /i(a>), and /2(c)) through a least square fit based on the steady state 
test data. 

21. The method of claim 13 wherein the step of determining the control yaw 
moment includes calculating the control yaw moment based on a derivative of the desired yaw 
rate. 

22. The method of claim 13 wherein the step of determining the control yaw 
moment includes calculating the control yaw moment based on a saturation function. 

23. The method of claim 13 wherein the step of determining the control yaw 
moment includes calculating the control yaw moment (M z ) based on the following equation: 

M z = JJdes - [*[(C FL +C fR )a f cos S + (?j FL + lFR )F zF sin 8] - b(C RR + C RL )a R + 

(C * C FL " d * °FR )a F SiT1 5 ~ C{11 FL F ZF cos 5 + ^RL F zR > + d ^FR F zF cos 6 + ^RR F zR M + KvSAT^j^-) 
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24. A method of estimating the retarding torque of an electromagnetic retarder 
having a control element, said method comprising: 

modeling the estimated retarding torque of the retarder (T es t) based on the following 
quadratic function of a velocity (co) of the control element and a magnitude of a retarder 
excitation current (i) 

estimating coefficient functions /o(co), /i(co), and /2(G)) from steady state test data 
performed for various rotor speeds. 

25. The method of claim 24 wherein the coefficient functions /o(co), /i(co), and 
/2(a)) are defined by 

f s (a)) = a i0 +a n a) + a i2 a) 2 

and wherein the step of estimating the coefficient functions /o(co), /i(co), and fi{(^) includes 
estimating parameters ay through a least square fit based on the steady state test data. 

26. The method of claim 25 wherein the step of estimating the coefficient functions 
/o(co), /i(cd), and /2(G)) includes recalculating the coefficient functions /o(co), /2(G)), and /2(G)) 
for each rotor speed and based on the estimates of parameters c^j. 

27. The method of claim 26 wherein the step of estimating the coefficient functions 
/o(co), /i(co), and fi{(b) further includes re-estimating the parameters osj based on the 
recalculated coefficient functions /o(co), /i(co), and /2(a)) through a least square fit based on the 
steady state test data. 



33 



28. An electromagnetic retarder assembly comprising: 

an electromagnetic retarder having a stator with conductive windings and a rotor 
rotatable relative to said stator; 

a sensor sensing a rotational speed of the rotor and generating signals indicative of the 
rotational speed; 

a controller communicating with said sensor to receive said signals, said controller 
configured to communicate an excitation current to said electromagnetic retarder and to 
estimate the retarding torque of the retarder (T es t) based on a quadratic function of the velocity 
of the rotor and the magnitude of the excitation current. 

29. The electromagnetic retarder assembly of claim 28 wherein the quadratic 
function is: 

T e5t =Mo)) + Ma))*i + f 2 (a))*i 2 

30. The electromagnetic retarder assembly of claim 29 wherein the coefficient 
functions /o(co), /i(a>), and / 2 (o)) are estimated from steady state test data performed for 
various rotor speeds. 

31. The electromagnetic retarder assembly of claim 30 wherein the coefficient 
functions /o(«), /i(co), and /2(co) are defined by 

f i (a>) = a i0 +a n a> + a S2 e> 2 
and wherein the parameters o&j are estimated through a least square fit based on the steady state 
test data. 
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32. The electromagnetic retarder assembly of claim 31 wherein the coefficient 
functions /o(co), /i(co), and / 2 (co) are estimated by recalculating the coefficient functions /o(o>), 
/2(a)), and /2(G)) for each rotor speed based on the estimates of parameters a X y 

33. The electromagnetic retarder assembly of claim 32 wherein the parameters a$j 
are estimated based on the recalculated coefficient functions /o(co), /i(co), and /2(a) through a 
least square fit based on the steady state test data. 

34. The electromagnetic retarder assembly of claim 29 wherein the controller is an 
open loop controller providing a current optimal torque without the need for a current feedback 
signal from the retarder. 

35. The electromagnetic retarder assembly of claim 34 wherein the controller is a 
parametric model control unit. 
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